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We explore time-symmetric hypersurfaces containing apparent horizons of black objects in a 5d 
spacetime with one coordinate compactihed on a circle. We hnd a phase transition within the 
family of such hypersurfaces: the horizon has different topology for different parameters. The 
topology varies from S'® to S^ x S®. This phase transition is discontinuous - the topology of the 
horizon changes abruptly. We explore the behavior around the critical point and present a possible 
phase diagram. 


Several black objects solutions exist in spacetimes with 
dimensionality greater then four containing compact di¬ 
mensions. Among these solutions are black strings (BS) 
and black holes (BH)- two black objects with distinct 
horizon topology. However no general analytic solution 
is known for a black object in a compactihed space-time 
in more then 4d spacetime. 

Gregory and Laflamme (GL) [1,2] discovered that a 
uniform black string, a product of a Schwarzschild solu¬ 
tion with a line, develops a dynamical instability if a com- 
pactihcation radius is ’too large’. They postulated the 
existence of a new branch of non-uniform black string so¬ 
lutions. The endpoint of an unstable uniform string was 
expected to be a BH. However, Horowitz and Maeda [3,4] 
argued that the GL instability cannot cause a uniform BS 
to decay to a BH, at least in a hnite affine time. This 
is because of the ’no tear’ property of horizons. They 
conjectured that the endstate of such a decay would be 
rather a non-uniform BS. Non-uniform BS solutions con¬ 
nected to the GL point were constructed perturbativly by 
Gubser [5]. Wiseman [6] found non-uniform BS solutions 
using fully non-linear numerical calculations. However, 
Wiseman stresses that these solutions are most likely un¬ 
stable. 

We consider here the hve dimensional spacetime with 
one compact dimension of length 2L. We denote the co¬ 
ordinate along the compact dimension by z. A neutral 
black object has a characteristic dimensionsfull parame¬ 
ter, its asymptotic 4d mass, M. It is convenient to de¬ 
scribe the system by a single dimensionless parameter, 
C = GiM/{2L), where the effective Newton constant is 
G 4 = G5/2L. If ^ «C 1 the black object is expected to re¬ 
semble a 5d Schwarzschild BH since local measurements 
cannot probe the compactness of the z direction. The 
horizon of this BH is therefore expected to have an 5® 
topology. At the opposite end, C ^ one expects that 
the horizon would extend over the compact dimensions 
wrapping it completely. The horizon’s topology of this 
object, S'^ X S'®, suggests the name ’BS ’. An uniform BS 
is described by a 4d Schwarzschild solution times a cir¬ 
cle. A non-uniform BS is characterized by a non-trivial 
z-dependence along the circle. An attempt, based on a 
thermodynamical reasoning, to construct the phase dia¬ 
gram that would include consistently different phases of 
solutions was recently undertaken by Kol [7]. 

In this letter we make a step toward the self-consistent 


study of black objects in compact 5d background. We 
solve numerically the initial-value problem for a moment 
of time-symmetry for spacetimes containing black ob¬ 
jects. We determine the apparent horizon of the existing 
black objects in these solutions and we follow the tran¬ 
sition from spherical to cylindrical topology in a family 
of solutions with a different dimensionless parameter C,- 
Although the conhgurations we hnd are not 5d static 
solutions, these solutions, give us some insight on the be¬ 
havior of the black objects in a compactihed background. 
A similar approach was previously applied [8] to explore 
the time-symmetric BH problem in the brane world sce¬ 
nario. 

We consider a time-symmetric slice through the black 
object’s space-time. To generate the black object solu¬ 
tion we consider a conhguration with an artihcial matter 
around the origin. This matter gives rise to a non-hat 
spatial metric along the slice. We solve numerically for 
this metric. Then we solve for the apparent horizon of 
the black object. The method allows us to determine the 
spatial geometry of the black object. Since we solve along 
a time-symmetric slice, this spatial geometry constitutes 
the initial data for subsequent dynamical time evolution 
of this black object’s geometry. Limiting ourselves only 
to the determination of the initial data we hnd explicitly 
time-symmetric solutions with a distinct horizon topol¬ 
ogy and construct the phase diagram for a family of such 
solutions parametrized by C- We show that there is a 
phase transition within this family. We identify the criti¬ 
cal (jo and show explicitly that if is reached from below 
a corresponding BH solutions becomes deformed and the 
horizon has a cigar-like shape. If <jc is approached from 
above, the corresponding BS solution becomes more and 
more non-uniform. However, this cannot be considered 
as evidence for existence of a non-uniform static BS solu¬ 
tions. This initial data may relax to a uniform string dur¬ 
ing its dynamical evolution. We hnd that the transition 
between both topologies is not smooth: the spherical-like 
horizon jumps suddenly to become a cylindrical-like hori¬ 
zon. Put differently, neither the BS pinches nor the BH 
has its south and north poles touching each other. 

The method that we describe here enables us to study 
properties of black objects with different apparent hori¬ 
zon’s topology in a single consistent numerical scheme. 
Our method does not require a prescription of the topol¬ 
ogy of the horizon. The topology is not predetermined 


1 


but rather it is a derived result of numerics. Usually, to 
obtain a static black object solution of the 5d Einstein 
equations one would had to specify the topology of the 
horizon. In other words there is no a single numerical 
scheme that can find black objects with distinct horizon 
topology. 

Let us consider a time-symmetric, spacelike hypersur¬ 
face T,t with a vanishing extrinsic curvature. The Hamil¬ 
tonian constraint reads 

, (1) 

where t^ is the unit normal to Sj and is the 5d stress- 
energy tensor. From now on we work in units where 
c = Gs = 1. The momentum constraint is trivially satis¬ 
fied provided that the matter is static and the extrinsic 
curvature of this slice vanishes. 

We choose the metric on Sj to be conformally flat 

dl^ = + dz^ + r'^dVL\) . ( 2 ) 

This ansatz has been adopted for simplicity. We didn’t 
expect that the static solution will be conformally flat. 
However, we expect that the trend we discuss is insensi¬ 
tive to this assumption. The method can be easily gener¬ 
alized to non-conformal choices. More complicated met¬ 
rics will be considered elsewhere. 

With this choice of the metric the constraint equation 
takes the form 

+ Pi’’^ = 0 , where = ^dr {^^dr) + dzz ■ (3) 

Here we have defined p = > 0, the physical en¬ 

ergy density. This p is the matter density, distributed 
around the origin. In practice, we pick a localized matter 
distribution that occupies some finite region around the 
origin. This artificial matter gives rise to a non-trivial ip 
in a way that does not involve specifying an inner bound¬ 
ary of a black object. For sufficiently concentrated mass 
an apparent horizon appears. For lower values of p the 
solution describes a momentary static star. 

To solve the elliptic Eqn (3) we have to specify bound¬ 
ary conditions. These conditions are: (i) At the equato¬ 
rial plane, 2 : = 0 , we have a reflection b.c. dzip = 0 ; (ii) 
At the symmetry axis, r = 0, we have also a reflection 
symmetry: drip = 0; (hi) Periodicity in 2 : implies also 
reflection symmetry at z = L: dzip = 0; (iv) At r —> 00 
we set Ip ^ 1 + G^M/r, which is consistent with the 
4d asymptotic behavior*. This condition is implemented 


*The asymptotic metric is expected to take the form of a 4d 
Schwarzschild times a line. This metric generally cannot be 
bronght to the conformally flat form (2). However, asymp¬ 
totically, as r ^ 00 , a transformation is possible at the lead¬ 
ing order. The asymptotic behavior of the conformal factor 
would become ip I + GaM/ r. The length of the circle tends 
asymptotically to 2L so we can read the asymptotic 4d mass 
in the units of the effective Newton constant G 4 = Gil2L. 


by rewriting it as dr (rip) = 1 , eliminating the need to 
specify M. Instead, M is determined from the numerical 
solution. 

There is a wide literature about equations like (3), see 

[9] and references therein. This equation is ill posed and 
it does not have a unique solution. In order to bring 
the equation to a well posed form we rescale the matter 
density p pip~^~‘^, with some s > 0 and nonphysical 
density p, see e.g. [9]. The resulting equation 

VV -b pV’”® = 0 (4) 

is well posed and has a unique solution. The rescalling 
of the matter density is unimportant as we are interested 
only in the external vacuum part. The solution for ip is 
found using relaxation. 

We obtain a sequence of momentary time-symmetric 
solutions by fixing the density of the artificial matter and 
its location as p = 1O®0(O.5 — r)0(O.5 — z), and varying 
continuously the length of the compact asymptotic circle, 
2L. The figures below are obtained for this source. We 
checked that our results are not affected by the specific 
choice of the source. Taking the smooth distribution p = 
10 ®exp(-rVo-^) [exp(-zVo-D + exp(-(z - 2 L)'^/al)] 
with various Ur and Gz we found the same overall pic¬ 
ture. The values of Cc varied, as expected, depending on 
the source. The variations in Cc are Cc = -98 — 1.8. More¬ 
over, taking Gr — 0.5, Gz 3> A, i.e. practically cylindrical 
source, we were able to reproduce the uniform BS solu¬ 
tion. The position of the horizon in this solution (see Eq. 
(5) and the subsequent discussion) was determined to 
within 0.1%. This provides an independed check on the 
overall accuracy of our calculation in addition to other 
standard tests of convergence, errors scaling etc. 

To solve the equation for ip we used grids, covering the 
domain 0 < z < L, and 0 < r < Rent, with typical grid 
spacings of Ar ~ 0.02 and Az ~ 0.01. The value of Rent, 
where the grid was cutoff and the asymptotic b.c. (iv) 
was implemented, has been taken as Rent = 5,10 and 20. 
We checked that the results are insensitive to variation 
of Rent, provided that Rent > 5. 

To envisage the spatial metric around the black objects 
in Fig.l we plot the contours of ip in two cases that corre¬ 
spond to a BH solution and a non-uniform BS solution. 
The matter is located near the origin and is encircled 
by the horizon in either case. The geometry outside the 
apparent horizon has an axisymmetric structure and it 
becomes asymptotically flat. The ip contour lines are 
spherical near the origin and become cylindrical as r in¬ 
creases. 

Once we obtain ip we determine the existence of the 
apparent horizon. An apparent horizon is defined by a 
zero-expansion of the null rays generating the horizon 

[10] . For the time-symmetric hypersurface Ej this condi¬ 
tion can be written as 

= 0 , (5) 

where is the the unit normal to the apparent horizon. 


2 



¥ 



r 

FIG. 1. Contours of ip. Upper panel: C ~ 1.85 and there is 
a non-uniform BS. Bottom panel: ~ 1.7 and there is a BH. 
In both plots one observes how the deformation of the contour 
lines fades asymptotically. We plot also the apparent horizons 
in both cases. These are designated by thick curves. In the 
upper plot that corresponds to the BS phase there are two 
horizons. The inner spherical apparent horizon, designated 
by the dotted thick curve, and the outer cylindrical horizon, 
designated by the solid curve. 

To simplify the treatment, we distinguish between two 
different topologies for the horizon. 

(1) When the horizon has the topology x we 
choose cylindrical coordinates (r, z) and we solve for a 
curve r = h(z). 

(2) When the horizon has the topology of it is conve¬ 
nient to transform to spherical coordinates R, x defined 
by r = i?sin(x), 2 ; = Rcos{x) The horizon in the R,x 
plane is given by a curve R = h(x)- 

The unit normal to the curve that defines the horizon 
is 


= (C, -Ch') . (6) 

The parameter C (r, z) could be read from the normaliza¬ 
tion condition = 1. In both cases we solve numeri¬ 
cally equation (5) to obtain the position of the apparent 
horizon. 

A useful qualitative parameter employed as a measure 
of the non-uniformity of a BS [5] is A = l/2(rmaa;/'Cmm — 
1 ) where and r^ax are the minimal and the max¬ 
imal 4d Schwarzschild radii of the apparent horizon. 
For a uniform string A = 0. In the BH phase A = 
oo. Therefore, for a BH we define another parameter 
A = Rrfiax /Rmin 1; where Rmax and Rjnin are the 5d 
Schwarzschild radii of the horizon. This parameter gives 
an idea of deformation of the BH’s horizon. 
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FIG. 2. The path in the configuration space of the numer¬ 
ical solutions. The BH (BS) phase is displayed in the bottom 
(top) panel. In both panels the horizontal dashed line that 
is located at — 1.78 designates the BS BH transition 
point. Above this line there is BS phase, while below the line 
there is the BH phase. The BS becomes uniform at ~ 3.0 
and the BH becomes spherical at ~ 0.15. 

We find that there are two topologically distinct ap¬ 
parent horizon solutions. At small C the topology of the 
horizon is and the horizon is close to be exactly spher¬ 
ical. When ( increases we see that the horizon begins to 
deform, deviating from a spherical shape but still remain¬ 
ing topologically 3-sphere. At a certain value, Cm — 1-78, 
a phase transition takes place - the topology of the ap¬ 
parent horizon changes form to x S^. In fact there 
are two apparent horizons in the BS phase. The outer 
horizon, has a cylindrical topology, while the inner one 
has a spherical topology. In Fig.2 we plot the sequence 
of solutions parametrized by C- In this figure there are 
finite values of C when A = 0 or A' = 0. In fact for 
these C the deformation of the horizon becomes so small 
that cannot be resolved by our numerics and we put the 
corresponding lambdas to zero. 

Another interesting result is the measure of geomet¬ 
rical deformation of the horizons. In the BS and the 
BH phases this measure is supplied by A and X' respec¬ 
tively. The non-uniformity of the BS is displayed in the 
upper panel of Fig. 2. Near the critical point the most 
non-uniform string has A cx 0.22. The non-uniformity 
disappear as ( increases. For C > 3.0 the BS becomes a 
uniform BS. The deformation of the horizon in the BH 
phase could be seen in the bottom panel of the same Fig¬ 
ure. The maximal radius, Rmax, always occurs at the 
axis, r = 0. The BH becomes more and more oblate and 
stretched along the symmetry axis, as we approach the 
critical point. The most deformed BH has A' ~ 0.15 just 
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before the transition. 

The phase transition is discontinuous. At the critical 
value of ('c the spherical horizon jumps suddenly to be¬ 
come a cylindrical. To get insight on the behavior of the 
phase transition, in the BH phase we have computed the 
proper distance along the r = 0 axis from the the BH 
horizon at the axis to z = L\ 

£ = [ '(/’(r = 0, zjdz . (7) 

J 

This distance decreases as we increase C,. One could ex¬ 
pect that as horizon grows and as £ ^ 0 the north pole 
of the BH would tend towards its south pole and they 
will touch. However, we find that as C —> Cc ^ reaches a 
finite value. We plot the behavior of £ as a function of 
C — Cc in Fig. 3. One observes that £ tends to a positive 
constant just before the transition. 
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FIG. 3. The proper distance from the BH to the reflection 
plane z = L does not decrease to zero but tends to a finite 
value as we approach the transition point. 

The initial data that we have constructed here is anal¬ 
ogous to the Misner initial data [11] for a family of mo¬ 
mentary static two equal mass BHs in 4d GR. Misner 
choose a sequence of conformally flat metrics with the 
conformal factor parametrized by a certain parameter /r 
that is related to the mass of the BHs and their proper 
mutual separation. As fj, varies the shape of the initial 
apparent horizons varies. If the BHs are close enough, 
that occurs for small /r, a new apparent horizon suddenly 
appears, surrounding both BHs on the initial hypersur¬ 
face. In other words, there is a critical /tq that divides 
two distinct possibilities for the topology of the apparent 
horizon on the initial slice. Just by looking at this ini¬ 
tial data sequence one has an indication that the event 
horizons for two BHs will merge and form a distorted 
BH during an actual evolution. The value of /i when 
this merge occurs, generally would not coincide with the 
theoretical, /tq. In fact, the numerical evolution [12,13] 
of Misner initial data shows that the qualitative picture 
obtained for the sequence of the initial data is correct. 
The actual critical value of /i does not coincide with /tq, 
they are, however, not that different from each other. 

Here we have an infinite array of BHs that are ap¬ 
proaching each other simultaneously. We have shown 
that there is a family of initial data parametrized by (^. 


When the separate BHs in the array are getting closer 
they become distorted form the spherical shape. At a 
critical value Cc the separate horizons are engulfed sud¬ 
denly by a single cylindrical-like horizon. The effective 
cylindrical horizon after the transition is non-uniform. 

It is important to stress that the sudden jumps of the 
apparent horizon topology are different from the first or¬ 
der transition that can happen between different phases 
of static solutions, discussed in [5,7]. This is because the 
apparent horizon that we discuss here is not casual as 
it is defined only locally. The event horizon is a global 
concept and it is expected to be larger than the apparent 
horizon. Only for static solutions both horizons coincide. 
Since our solution is not static the jumps of the appar¬ 
ent horizon cannot exclude the possibility of a smooth 
transition between the BHs and the BSs event horizons, 
discussed in [14]. 

The concrete numerical value of Cc isn’t important as 
it is just a number characteristic to the specific initial 
data sequence. However we expect that the qualitative 
behavior would be similar in the static solutions as well. 
We believe that a dynamical evolution of our initial data 
would confirm this qualitative picture and will yield ac¬ 
tual critical value, (c- 
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